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Abstract — In the control of classical mechanical systems, feed- 
back has been applied to the generation of desired nonlinear 
dynamics, e.g., in chaos control. However, how much this can 
be done is still an open problem in quantum mechanical 
systems. This paper presents a scheme of enhancing nonlinear 
quantum effects via the recently developed coherent feedback 
techniques, which can be shown to outperform the measurement- 
based quantum feedback scheme that can only generate pseudo- 
nonlinear quantum effects. Apart from the advantages of our 
method, an unsolved problem is that the decoherence rate is also 
increased by the quantum amplifier, which may be solved by 
introducing, e.g., an integral device or an nonlinear quantum 
amplifier. Such a proposal is demonstrated via two application 
examples in quantum optics on chip. In the first example, we show 
that nonlinear Kerr effect can be generated and amplified to be 
comparable with the linear effect in a transmission line resonator 
(TLR). In the second example, we show that by tuning the gains 
of the quantum amplifiers in a TLR coherent feedback network, 
the resulting nonlinear effects can generate and manipulate non- 
Gaussian "light" (microwave field) which exhibits fully quantum 
sub-Poisson photoncount statistics and photon antibunching phe- 
nomenon. The scheme opens up broad applications in engineering 
nonlinear quantum optics on chip. Particularly, in this study, 
the concept of feedback nonlinearization which is very useful 
for quantum feedback control systems is introduced. This is in 
contrast to the feedback linearization concept used in classical 
nonlinear feedback control systems. 

Index Terms — Feedback nonlinearization, quantum coherent 
feedback control, nonlinear quantum optics, on-chip quantum 
optics, quantum control. 



L Introduction 

OVer the last decades, the control of quantum phenom- 
ena has been steadily advanced in many fields such 
as quantum communication and computation, laser-induced 
chemical reaction, and nano electronics [1|-|39|. However, 
the implementation of realtime feedback, which is the core 
of control theory and engineering, is still in its infancy. The 
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major technical obstacles are (1) the time scale of general 
quantum dynamics is too fast to be manipulated in realtime 
by currently available electronic devices; (2) the required 
quantum measurements are generally extremely hard to do; 
and (3), more essentially, the back action brought by the 
quantum measurement deeply annoys the control designers 
because it keeps dumping entropy into the system before the 
feedback attempts to reduce it. So far, it is still not clear to 
what extent quantum feedback control may outperform the 
open-loop control, and this impedes the discoveries of new 
applications of quantum feedback control. 

Up to now, there have been two commonly studied classes 
of quantum feedback strategies in the literature. The first one 
is called the measurement-based feedback control |401-||46t. 
A typical implementation of such a scheme is to shoot an 
electromagnetic probe field through the quantum system to be 
controlled, carrying a part of information of the system that 
can be detected by some measurement apparatus and converted 
into classical signals, which are then fed back to adjust the 
input of the system. This is the analog of classical feedback 
loop. However, the accompanied problem is that the quantum 
measurement disturbs the system (i.e., the measurement back- 
action), and thus adds unremovable noises. Such a feedback 
was shown in our study BTl to be only capable of generating 
"classical" nonlinearity, which is not fully quantum. 

By contrast, the other strategy is called coherent feedback 
control [i48l - ll52l . in which the probe field is, instead of being 
read out, coherently guided back into the system after being 
unitarily transformed via quantum controllers (e.g., quantum 
beam splitters, quantum switchers, and quantum amplifiers). 
Such a strategy preserves the quantum coherence of the 
system, and is completely new to the theory of control system. 

This paper will propose a promising application of quantum 
feedback control that is capable of inducing fully quantum 
nonlinear effects into the controlled systems. Similar idea has 
been applied in the pioneer work |[53l for a novel approach to 
engineer the nonlinear dynamics of the nanomechanical sys- 
tem by nonlinear feedback control. This is important because 
nonlinear quantum processes [54] are essential in engineering 
many interesting quantum phenomena, such as the photon 
blockade induced by nonlinear Kerr effect and the generation 
of non-Gaussian light, with broad applications in quantum 
information processing, quantum nondemolition measurement, 
and the preparation of particular nonclassical states, e.g., the 
Schrodinger cat state. However, the nonlinear effects induced 
by the natural field-matter interactions are normally very weak. 
Therefore, artificial enhancement of quantum nonlinearity is 
crucial. Note that a similar idea has been employed in classical 
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control of chaos, but the case for quantum control has been 
rarely studied f53l, and, as will be seen below, is much more 
complicated due to the so-called quantum coherence. 

As a potential application, our coherent feedback nonlin- 
earization scheme can be applied to the emerging field of 
quantum optics on chip, i.e., producing optical-like quantum 
phenomena on compact solid-state chips such as waveguide 
circuits or superconducting circuits fSSj-fSSl. The on-chip 
optical devices are powerful in demonstrating particular optical 
phenomena that are hard to be observed in conventional optical 
setups. So far, the existing on-chip optical experiments are 
mainly done in the linear regime, due to the fact that the 
natural nonlinear effects induced by the couplings between the 
optical fields (the microwave fields) and the solid-state devices 
are too weak to be observed. In this paper, we will show that 
such nonlinear effects of on-chip lights can be "artificially" 
generated and amplified via the coherent feedback strategy. 

This paper is organized as follows. In Sec. [Ill preliminaries 
are given for a brief introduction of the theory of the co- 
herent feedback control network. In Sec. |III] we introduce 
the basic setup and the mathematical model of the coherent 
amplification-feedback loop, from which the Hamiltonian of 
the controlled system can be effectively reconstructed. In 
Sec. IIVI we apply these general results to the generation of the 
Kerr effect and the cross Kerr effect in the superconducting 
circuits. Furthermore, in Sec. |V] we study how to construct 
general fourth-order controllable nonlinear quantum Hamilto- 
nians and their applications to nonclassical on-chip "lights" 
(microwave fields). Conclusions and perspectives are given in 
SeclVj 

II. Preliminaries 

In quantum mechanics, the state of an isolated system can 
be described by a vector \ip) in an abstract Hilbert space 
H, on which the system observables can be described by 
operators. The evolution of the system state is governed by 
the Schrodinger equation: 

\i;{t))^-iH\^{t)), 

where the Hamiltonian H is an Hermitian operator denoting 
the "energy" observable of the quantum system. Here, the 
Planck constant h has been assigned to be 1. If the system 
is bathed with the external environment, then, instead of the 
state vector lip), the system state should be represented by the 
so-called density operator p, which is Hermitian and positive 
semi-definite on H. The evolution of p can be described by 
the following master equation [|59l : 

p{t) = -i[H,p{t)] + CD [p{T)\Te[o,t]], 

where the commutator [•, •] is defined as [A, B] — AB — BA. 
The superoperator Cd [p (t) \t E [0, t] ] represents the dissipa- 
tion in the system due to the interaction with its environment, 
which depends on the system states in the whole time interval 
[0,t]. Under the so-called Markovian approximation to omit 
the back action effects from the environment, we can obtain 
a time-local dissipation superoperator Cd [p {t)]- 
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Fig. 1. (Color online). The schematic diagram of the input-output system. 



Generally, in quantum optics, the quantum electromagnetic 
field (e.g., the probe light to be used below) is treated as the 
collection of the quantized modes whose Hamiltonian reads 



Hp. = 



where is the annihilation operator of each quantized field 
mode with frequency w. Each b^ acts on the corresponding 
Hilbert state space spanned by the Fock states \n^), riuj = 
0, 1,2, - -, where \nuj) represents the state of the quantized 
field mode that contains ricj photons with frequency oj, in the 
way that buj\ni^) is proportional to \nuj — 1). In addition, b^j 
also satisfies the continuous-variable canonical commutation 
relationship b^, bl^, = S {uj — uj'). Traditionally, we use the 



time-dependent operator 

m = 



1 



diue 



to represent the field consisting of a continuum of field modes, 
which satisfies [b{t),b'' {t')] =6{t-t'). 

Next, we give a brief review of the theory of the quantum 
feedback network developed recently fST], ll52l . Il60l - ll70l . 
Consider the general input-output system given in Fig. [T] 
The input of the quantum plant (with internal Hamiltonian 
H) bin(t) — [bi{t), ■ ■ ■ , bn{t)]'^ contains n mutually different 
fields bi{t),i = 1, all initially in vacuum states. 

The input field bin (t) transmits through a quantum beam 
splitter described by an n x n unitary scattering matrix S, and 
then interacts with the plant through the dissipation channels 
represented by the dissipation operators L = [Li, ■ ■ ■ ,Ln]^. 
In the following discussions, we will concentrate on the single 
input case, i.e., n = 1, to simplify our discussions. 

Let V{t) be the unitary evolution operator of the composite 
system of the plant plus the input field (as a quantum system). 
Correspondingly, the control system can be described by the 
following quantum stochastic differential equation ||68l . ITTTI : 



dVjt) 
dt 



bl{t)iS-I)Vm^{t) + bl{t)LVit) 

1 



^L^SV{t)binit) 



-L'^L-iH 



V{t), (1) 



with initial value Vb = /, where / is the identity operator The 
output field is [72\: 



boutit) = V{t)^SV{t)hn{t) + V{t)^LV{t). 



(2) 



Equation ^ is the Wick-ordered differential equation of Vt 
(creators appear on the left, annihilators on the right), which 
is equivalent to its quantum stochastic differential form (see, 
e.g., Eq. (30) in Ref. |I52)- Here, L and H (H is self-adjoint) 



SUBMITTED TO SPECIAL ISSUE OF IEEE TRANSACTIONS ON AUTOMATIC CONTROL 



3 



(Sj, Lj, Hj) (S„ L„ H,) 



K 




o- 




















o— 





Fig. 2. (Color online). The schematic diagram of the series product system. 
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Fig. 3. (Color online). The schematic diagram of the coherent feedback 
system. 

are system operators commuting with fields (s) and bj^^ (s) 
for earlier times s < t, which means that they are adaptive. 
S is unitary satisfying S^'S ^ 55^ = /, which, for the one- 
dimensional case we consider {n — 1), is just a phase factor 
S — exp(z0). In a compact notation, the above input-output 
system can be represented by {S, L, H). 

After averaging over the vacuum input field 6in(i) which 
can be treated as a quantum random process, the system given 
by Eq. ([T]i can be transformed to the master equation for the 
density operator of the plant [52 1: 

p{t) = ^i[H,p{t)]+V[L]p{t), (3) 

where the superoperator T) [L] p is defined as: 

V[L]p = LpL^ -]^L'^Lp-]^pL^L. 

Let us consider the two cascade systems shown in Fig. |2] 
By introducing the Markovian approximation to omit the time 
delay for the output of the first system (5*1, Li,Hi) to reach 
the second system (5*2, L2, H2), the total system in the series 
product can be described as follows |52|: 

(^S2Si,L2 + S2LuHi+H2 + ^(^L\sIL2-LIS2Li^^ . 

(4) 

As a special case, if we feed the output of the system 
{S,L,H) back and take it as the input of the same system 
to construct a direct coherent feedback network as shown in 
Fig. [3] From Eq. such a feedback network can be described 
by: 

(^S^,L + SL,H+^L^S'' ~S)L^ . (5) 

III. Quantum amplification-feedback loop 

In order to fulfil strong and controllable feedback-induced 
nonlinear effects, we consider the modified quantum coherent 
amplification-feedback loop shown in Fig. |4] There are two 
differences between the traditional coherent feedback network 




Plant 

Fig. 4. (Color online). The schematic diagram of the coherent feedback and 
amplification network. 

given in Fig. [3] and the quantum coherent amplification- 
feedback loop given in Fig. |4] Firstly, the dissipation operator 
L which represents the interaction between the system and the 
input field 6in may be different from L / which represents the 
interaction between the system and the feedback field &out- 
Secondly, we add a quantum amplifier f73\ in the feedback 
loop, and feed the output field bout{t) from the system into it. 

The simplest setup to fulfil such a quantum amplifier is 
a driven squeezed cavity field [74] with the damping rate 
K and the tunable squeezing coefficient ^, which can be 
realized by the strategy given in Refs. flS], f76\ (see Fig. |5]l. 
As presented in Refs. Il75l . [|76[, the cavity field is coupled 
with a three-level atom which is further driven by a classical 
field. By adiabatically eliminating the degrees of freedom of 
the auxiliary three-level atom, we can obtain a controllable 
squeezed field in the cavity in which the squeezed coefficient 
is tunable by changing the coupling strength between the 
classical driving field and the three-level atom. This strategy 
can be extended to the solid-state superconducting circuit by 
replacing the cavity by a transmission line resonator (TLR) and 
the three-level atom by an auxiliary flux qubit (see Sec. V in 
Ref. ITTTI ). which is used to construct an on-chip amplifier in 
the following discussions. 

The Hamiltonian of the controllable squeezed cavity field 
can be represented, in the rotating frame, as: 

ffc = I (c^ 2 - c2) + {e't'c + c^e-'t') , (6) 

where c is the annihilation operator of the cavity field; 
A, 4> £ TZ represent the normalized amplitude and the initial 
phase of the classical control field driving the cavity mode. 
In the {S,L,H) notation, the squeezed cavity system can be 
represented by: 

(^I, ^c, I (c^' - c') + ^A {e'^c + c^e-'^^)^ . (7) 

The original system (5*, L, H) and the squeezed field can be 
looked as a series product system, which can be represented 
from Eq. (|4]i as: 

{s,^/^c + L,H + H,+ ^-^/l^{L'^c^c^L)^. (8) 
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Fig. 5. (Color online). The single-mode squeezed field which works as the 
quantum amplifier. 



The total coherent amplification-feedback loop given in Fig.|4] 
can be looked as a series product system of the subsys- 
tem (5*, L, H), the quantum amplifier, and the subsystem 
{S,Lf, H), and thus can be described by: 



{S^,Lf + S {y^c + L),H + H, + ^^{L^c - c^L) 
(it + ^ct) S^Lf - l\S {L + v^c) 



(9) 



Following the ideas of Refs. Il78ll . Il79l . we can adiabatically 
eliminate the degrees of freedom of the cavity mode by the 
singular perturbation approach to obtain the following master 
equation (see the derivations in Appendix): 



dp = 



H + cosh (ro) ( ^l'^jSL - ^L'' S'' L 



sinh(ro) ( (l^ - l\s) (i^ + l\S 



--Ae^'f' [(cosh(ro) + l)(L + 5ti/) 
+ sinh(ro)(it+it5^j +h.c.}, 
+Vs [L- S''Lf]p + V[L- S^'Lf] p, 



(10) 



where 

v.. 



L 



P = 



L 



NV 



{N + 1)V 

+M* [LpL - Dp/2 - pD/2 
+M (DpD - L^^p/2- pD^/2 



and 

ro = In 



N 



cosh (2ro) — 1 



M 



sinh (2ro) 



The two dissipation channels 2?[L — 5^L/]p and 
2?s [i — p are induced by the vacuum field bin 

and the squeezed vacuum field 

= cosh (ro) Shn ~ sinh (ro) &^„S'^ 

where 6?^ is generated by the quantum amplifier from bin, 
which satisfies that 

btnit)bt!it') - iN + l)6it~t'), 

btni^btnif) - M*6{t^t'), 
bt^it)bt!it') = MSit-t'), 



Under the adiabatic approximation, the input-out relation of 
the squeezed component can be written as: 



^out = vGo&out + \/Go — lb, 



t 

'out • 



It can be seen that the squeezed component in this case can 
be looked as a phase-insensitive quantum amplifier ||73l , ||74l 
with power gain: 



Go = cosh^ (ro) 



e 



(11) 



The system we discuss cannot be expressed by the {S, L, H) 
notation given in Ref. [52 1 due to the existence of the squeezed 
bath term Vg \L — S^Lf \ p. However, it can be checked that 
Eq. ( [Tol l coincides with those equations in Ref. |69| for linear 
quantum systems, in which linear quantum dynamical network 
elements including static Bogoliubov components (such as 
squeezers lISOl . lIHTl discussed here) are formulated by the 
transfer function and input-output equation. 

If the power gain of the quantum amplifier Go is far greater 
than 1, the master equation ( fTOl i can be simplified as: 



[ffeff,p]+GoI? 



(12) 



where 



i?cff = 



+ ■v/Go ^ cos 



2 f 2 



^ ^L^ + L]S) +h.c. 
L + + S'^Lf + l\.s) (13) 



represents the effective Hamiltonian under coherent feedback. 
One can immediately see that the system Hamiltonian has been 
reconstructed by the feedback loop involving the following 
tuple of parameters: 

C = {S, L, Lf, Go, A, cf,}, 

which can be properly designed to realize desired quantum dy- 
namics in the closed-loop system. Note that these parameters 
can also be chosen to be time-variant to get more flexibility, 
but such a case will not be discussed in this paper 

A system is said to be linear if its Hamiltonian H as a 
polynomial of the annihilation and creation operators a and 

is up to the second-order and the dissipation operator L is 
a linear combination of a and a^, otherwise it is said to be 
nonlinear. Obviously, in Eq. (|9]l, nonlinear dynamics can be 
generated by 



Hnl — 



'Go 



Lis 



L 



h.c. 



Go ( 1-L\SL - l-L^S^L 



2 f 



(14) 



if X or Z// is a second-order or higher-order polynomial 
of the annihilation and creation operators. What we have 
discussed above is the concept of feedback noninearization. 
More importantly, strong nonlinear effects can be produced 
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provided that the power gain Go of the quantum ampUfier in 
the coherent feedback loop is sufficiently high. 

In practical experiments, the power gain of the quantum 
amplifier cannot be too large. For example, in optical systems, 
the quantum amplifiers are typically implemented using a 
nonlinear optical material. Thus, the production of a high 
gain amplifier may require an optical material with strong 
nonlinearity which is hard to be realized. One possible solution 
of this problem is to cascade a series of low gain quantum 
amplifiers to obtain a high gain just like what we have done for 
classical systems. There is also a great progress for nonlinear 
amplification in solid-state quantum systems such as super- 
conducting circuits. For example, as shown in Ref. ll82l . the 
authors found that quantum signals can be greatly amplified 
near the bifurcation point of a nonlinear quantum device, 
which has been improved by the succeeding experiments. 
Thus, it is possible to realize a quantum amplifier with high 
gain in solid state systems other than optical systems. 

The idea of the above feedback nonlinearization strategy 
can be demonstrated in the following simple model. Let 
the plant be a single-mode field with normalized position 
and momentum operators as x — (a + a^) /\/2 and p = 
{—ia + ia^) /\/2- The internal Hamiltonian of the field is 
H — wfltfl. We choose L — ^J^x^, — ™d S = e*'^/^ 
respectively. If we coherently feed back the quantum signal 
(i.e., being an operator): 

where bin is the input vacuum noise, then i7„/ given in Eq. ( fT4l l 
contains third-order nonlinear terms of x and p. 

In comparison, the feedback signal in the measurement- 
based quantum feedback scheme ||£7l is a classical signal (i.e., 
being a scalar): 

kutit) = v^{x^)+m, 

where (A) — trAp is the average over the system state p, and 
^{t) is a classical white noise satisfying: 

It has been shown in Ref. BtI that such a feedback loop 
can only introduce "classical" nonlinearity for the controlled 
quantum system. In fact, statistically, only the trajectories 
of the expectation values of the operators x and p behave 
nonlinearly. Those of the higher-order quadratures remain 
linearly, which can only be altered by fully quantum nonlinear 
Hamiltonian terms such as Hni introduced by the coherent 
feedback. Such fully quantum nonlinear dynamics is essential 
to important physical applications such as light squeezing and 
the generation of the Kerr effect. 

IV. Generation of Strong Kerr Effects 

In this section, we will introduce the coherent feedback 
to generate strong and controllable nonlinear effects in the 
following two systems: 
(1) Kerr effect in a single-mode field with the annihilation 
operator a, where the nonlinear Hamiltonian to be con- 
structed is ffxcrr — Xi^^c) T X ^T^'-> 




Fig. 6. (Color online). Schematic diagrams of the superconducting circuits 
to generate strong and controllable Kerr and cross-Kerr effects by coherent 
feedback and feedforward with (a) for Kerr nonlinearity and (b) for cross-Ken' 
nonlinearity. 

(2) Cross Kerr effect in a two-mode field with the annihila- 
tion operators a and b, where the nonlinear Hamiltonian 

is i/cross-Kcrr = Xab {a'^ o) {b'^ b) , Xab € H- 

The generation of strong Kerr and cross-Kerr effects is 
crucial to nonlinear quantum optical phenomena, and has 
important applications to the generation of particular quan- 
tum states, e.g., the Schrodinger cat state |83| and universal 
quantum computation |84|, i.e., the construction of two-qubit 
CNOT gate. Let us see how the Kerr effect can be generated 
in on-chip quantum optics ll55l -ll58l realized by the super- 
conducing circuit Il85l - ll87l shown in Fig. |3a). The vacuum 
input field bin transmits through a 7r/2 phase shifter which can 
be implemented by the on-chip quantum beam splitter pSS] 
proposed recently. Then, the output field of the phase shifter 
is coupled to the fundamental mode of the electric field in a 
TLR via a charge qubit. 

The Hamiltonian of the coupled qubit- TLR system can be 
expressed as |[85l - ll87l : 

HqT — 4:Ec {n — Ug)^ — 2Ej cos ( ) coscj) + Waft^a, 

(15) 

where cJq and a are the angular frequency and the annihilation 
operator of the electric field in the TLR; (p is a phase operator 
denoting the phase drop across the superconducting loop of 
the charge qubit; n — —id/dcj) is the conjugate operator of 
(f), which represents the number of Cooper pairs on the island 
electrode of the charge qubit; the reduced charge number Ug 
on the gate of the charge qubit, in units of the Cooper pairs, 
can be given by ng = —CgVg/2e; Cg and Vg are the gate 
capacitance and gate voltage; Ec = e^/2 {Cg + 2Cj) is the 
single-electron charging energy of the charge qubit; E'j and 
Cj represent the Josephson energy and the capacitance of a 
single Josephson junction; and $o is the quantum flux, in 
Eq. ( fTSl ) denotes the external flux piercing the SQUID loop of 
the charge qubit which can be expressed as: 

^x^<^e+ VT (a + at) + 77i„ (bin + foL) : (16) 

where $e is the flux generated by the classical magnetic field 
through the SQUID loop, and rjT, Vin have units of magnetic 
flux and their absolute values represent the strengths of the 
quantum flux in the SQUID loop induced by the electric 
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flux qubit 




Fig. 7. (Color online). Schematic diagram of the quantum amplifier in 
superconducting circuit. 



field in the TLR and the input field 6i„. Let Ug — 1/2, 
<i>e = — $0/2, and adiabatically eliminate the degrees of 
freedom of the charge qubit, we can obtain the following 
effective Hamiltonian to represent the coupling between the 
TLR and the input field (see the derivations in Appendix): 



: UJaO'a 



$3 



(17) 



Here, we have omitted the linear term of a and a^, which can 
be compensated by a classical driving field imposed on the 
TLR. 

With the {S,L,H) notations, such a system can be ex- 
pressed as: 

H = ujaa^'a, L = ^7^0+0, S = e"/^ 



where 



7q 



$6 



is the damping rate of the TLR induced by the input field. 
The output field of the system is fed into a quantum amplifier, 
and then fed back into the plant via the damping channel 
represented by 



As analyzed in Sec. [Till the quantum amplifier used here 
can be implemented by a squeezed cavity field with a tunable 
squeezed coefficient ^, which can be realized by the supercon- 
ducting circuit given in Fig. |7] In this circuit, a transmission 
line resonator is driven by the quantum input field fcout and a 
classical driving field 

/ {t) ~ At cos {ujTt + 4>t) , 

where the angular frequency of the driving field is equal to 
the angular frequency ut of the fundamental mode of the 
electric field in the quantum amplifier, i.e., the TLR. The 
output field of the quantum amplifier is 6out- The TLR is 
coupled with a flux qubit which works as a A-shaped three- 
level artificial atom |89|. As discussed in Sec. V of Ref. llTTl . 
by adiabatically eliminating the degrees of freedom of the flux 
qubit, the effective Hamiltonian of the quantum amplifier can 
be written under the rotating wave approximation as: 



HcS = (jJtc'c + ^ 



where c is the annihilation operator of the quantum amplifier, 
and ^, 51, and i/j are tunable parameters. Let 57 — 2ujt, 
ip = —tt/2, the effective Hamiltonian Hctf can be written 
in the interaction picture as He given in Eq. (|6]l. With the 
experimentally realizable parameters, the squeezed coefficient 
^ can be as large as the damping rate of the quantum amplifier 
K (see, e.g., Ref. ITTTI ). Thus, the quantum amplifier with large 
power gain Go given in Eq. (fTTT i can be obtained if we tune 
the squeezed coefficient ^ such that ^ w k. 

With this setup and under the condition that Go ^ 1, an 
nonlinear Hamiltonian of the closed-loop dynamics of the su- 
perconducting circuit shown in Fig. |6ia) can be reconstructed: 



(18) 



where 6 — 2ATvGo7a and x — ^.y/Gola are the angular 
frequency shift and the strength of the nonlinear Kerr effect 
induced by the coherent feedback control respectively. As 
shown in Eq. ( fTSl l, the Kerr effect is enhanced by increasing 
the power gain Gq, which can be done with the on-chip ampli- 
fication device f73l. As a numerical example, if the parameters 
Wa/STT = 500 MHz, 7a/27r = 1 MHz, A|/27r = 576 MHz, 
and the power gain of the quantum amplifier Go = 100, then 
it can be calculated that 

(wa - S) /27r = x/27r = 20MHz. 

The strength of the generated Kerr term, which is comparable 
with that of the lower-order term, is about 10''^ — 10^ stronger 
than the Kerr effect induced by the natural coupling between 
the electric field in TLR and the nonlinear element in su- 
perconducting circuit (only around tens of kHz). Thus, with 
experimentally realizable parameters, the coherent feedback 
strategy may dramatically enlarge the nonlinear Kerr effect. 

Furthermore, using the same idea, if the plant includes an- 
other superconducting circuit whose TLR annihilation operator 
is b (see Fig. |6lb)), we can feed the amplified output field into 
this circuit via Lf ~ y/^b'^b, and a cross-Kerr Hamiltonian is 
obtained as below: 



s-Kc 



Xaba 



'ab%, 



whose strength Xab = ^.^JGo^alb can also be enhanced by 
increasing the amplification gain Gq. 

V. Controllable Fourth-order Nonlinear 
Dynamics 

This section will focus on the design of a controllable 
fourth-order TLR Hamiltonian: 



Hc?L = UJaO^a + ^ Xkx'l, 
k=l 



(19) 



where x., 



+ a) / \f2 is the normalized position operator 
of the TLR; and Xfe, = 1,2,3,4 are the coefficients of the 
fc-th order quadratures, which are all tunable parameters. iJoff 
given in Eq. ( fT9] l can be used to produce more interesting 
nonlinear quantum effects. The terms in Eq. (fT9] l have dif- 
ferent applications, e.g., the X2-term can be used to realize 
controllable squeezing in TLR ||62l ; the xa-term can be used to 
construct the cubic phase gate which is fundamental to realize 
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Fig. 8. (Color online), (a) Schematic diagram used to generate the fourth- 
order controllable Hamiltonian given in Eq. l ll9t . (b) The A-shape transition 
of the auxiliary flux qubit. 



universal continuous variable quantum computation | |90| ; and 
the X4 term is useful for generating the Kerr effect. More 
importantly, as shown below, the nonlinear Hamiltonian given 
in Eq. ^ can be used to generate non-Gaussian "light" (mi- 
crowave field) to show fully quantum sub-Poisson photoncount 
statistics and photon antibunching phenomenon |54|. The non- 
Gaussian "light" generated is possible to be used to transmit 
quantum information, which may have higher capacity of the 
information transmission than the Gaussian light in continuous 
variable quantum communication 191]. 

Such a nonlinear Hamiltonian can be constructed via the 
superconducting circuit in Fig. |8ja), where a TLR is coupled 
with a flux qubit (left) and two charge qubits (central, and 
right). The central charge qubit works as an auxiliary device 
for the quantum detection. By tuning the parameters of this 
auxiliary charge qubit, we can execute a detection of the 
square of the normalized position operator Xa by the probe 
field through the TLR. With this detection, we can obtain a 
quadratic damping operator L ~ ^x\, where 7 is the related 
damping rate. This is similar to the quantum measurement 
strategy given in Ref. |92| to detect the square of the nor- 
malized position operator of a nanomechanical resonator The 
charge qubit on the right plays the same role as in Sec. |IV] 
to induce a nonlinear Kerr term. The flux qubit in the circuit 
can be treated as a three-level artificial atom with A-shape 
transition ll89l . Il93l . whose interaction with the electric field 
in the TLR induces another nonlinear Hamiltonian. As shown 
in Fig. [Itb), a two-photon exchange process occurs between 
the cavity mode a in the TLR and the travelling-wave mode 
62 via the flux qubit, i.e., the A-shaped three-level artificial 
atom. Two photons in the cavity mode a with the same angular 
frequency uja annihilate, and one photon in the travelling-wave 
mode 62 is created in this process, and vice versa. Finally, the 
output field is coherently fed back after amplification to drive 
the electric field in the TLR, which leads to a third nonlinear 
Hamiltonian. 

Mathematically, the above setup results in three feedback 
loops with 

^1 = 5*2 = ^3 - e'-/2^ 
^1 = ^2 = ^3 = yfl^'i^ 

The corresponding parameters Aj 



72 a 



t2 



L 



3/ 



inputs of the quantum amplifiers and the phase shifters are 
tunable parameters need to be designed. Here, we let (pi — 
(f>2 = 0, 03 = — TT, 04 = — 7r/2. From Eq. ( fT3] l, we can obtain 
the desired effective Hamiltonian (T% with 



X3 = 2VG3773, X4 = 2^61771, 

where we have set parameters G2 — 0171/72 and A2 — 
Ai-\/ 72/71 in order to obtain the Hamiltonian form shown 
in Eq. (fT9] l. Therefore, by tuning the control parameters 
Gi, G3, Ai, A3, and A4, we can independently change the 
coefficients Xk, k — 1,2,3,4. 

These nonlinear terms can be designed to generate nonclas- 
sical microwave field (i.e., the so-called nonclassical "light") 
in TLR. To illustrate the effectiveness of the coherent feedback 
scheme, we set the parameters: 

Wa/27r ^ 100 MHz, 7i/27r ^ 73/27r = 1 MHz, 
7/27r = 1 MHz, d = G3 = 10^, 
A2/27r = 40MHz, A§/27r = 152.1 MHz, 
Al/27r = 200 MHz. 

As shown in Fig. |9|a), we can observe the sub-Poisson 
photoncount statistics indicating by the fano factor 

F^{{N'j~{Naf)/{Na)<l, 
and the photon antibunching phenomenon indicated by 

.g(2)(r)>g(^) (0) , 

where Na = a^a is the photon number operator of the TLR; 
(•) is the average over the system state; and the normalized 
second-order correlation function g^'^^ (r) is defined by; 



(2). N ^ (gt jt) gt {t + T)a{t + T)a (t)) 



1,2,3,4 of the 



Here a (t) is the operator in the Heisenberg picture defined by 
tr(a(i)po) — tr(a/9(t)), where po and p{t) are the initial 
state of the system and the state of the system at time t. 
Sub-Poisson photoncount statistics and photon antibunching 
phenomenon are typical quantum phenomena which violate 
the Cauchy-Schwartz inequality for the classical lights. Ad- 
ditionally, different from the Gaussian lights, e.g., the laser, 
which are quite similar to classical lights, the on-chip light 
generated is indeed a non-Gaussian light, which is highly 
nonclassical. In Fig. |9tb), we use the measure 

^^p^ = — ^] — 

to evaluate the non-Gaussian degree of the light generated ll94ll . 
where cr is a Gaussian state with the same first and second- 
order quadratures of the non-Gaussian state p. Simulation 
results in Fig. |9jb) show that high-quality non-Gaussian state 
with S [p] > 0.25 can be obtained. As pointed out by Ref. ||95]| , 
the maximal value of S [p] is not larger than 1/2 for single- 
mode quantum states. Thus, the non-Gaussian degree of the 
generated light can be larger than half of the maximal value 
that can be reached by any non-Gaussian states. 
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(a) 



F<1 Sub-Poisson 
' statistics 



■ 1 2 3 4 5 6 




(c) 5^^^=0.27 




/ V[p]>0 non-Gau 


ssian light 



Fig. 9. (Color online) (a) Fano factor F{t) and (b) the normalized second- 
order correlation function ^'■^'{t): the sub-Poisson photoncount statistics 
indicating by F < 1 and the photon antibunching indicating by (;'^'(t) > 
g^'^") (0) seen in small time scale exhibit the nonclassical properties of the 
light generated, (c) non-Gaussian measure S [p]: a highest value 0.27 can 
be attained, and it is shown that the stationary state of the electric field is 
also non-Gaussian and thus nonclassical state, t* = 0.2 ns is a normalized 
time-scale. 



Finally, it should be pointed out that the decoherence 
(i.e., the destruction of quantum coherence) rate will also 
be increased by the quantum amplifier that is designed to 
enhance the nonlinear effect. Thus, there is a tradeoff between 
pursuing strong nonlinear effect and weak decoherence, which 
may somewhat limit the ability of effectively manipulating 
nonlinear quantum optical phenomena, e.g., the generation of 
Schrodinger cat state. This problem will hopefully be solved 
by introducing a nonlinear amplifier whose signal/noise ratio 
is high. 



VI. Conclusion 

In summary, we present a method of engineering strong 
and controllable nonlinear effects in quantum systems by 
coherent feedback control and amplification. A byproduct of 
this investigation is the introduction of the concept of feedback 
nonlinearization which is very useful for quantum feedback 
control systems. To the authors' knowledge this concept has 
never been discussed in the literature. The applications in TLR 
superconducting circuits demonstrate its power of generating 
strong and controllable nonlinear Kerr and cross-Kerr effects 
and more complex nonlinear phenomena. They open up new 
perspectives to the design of nonlinear circuits for quantum 
optics on chip, where the systematic design methodology of 
the feedback loop parameters, including S, Gq, etc., for more 
complex nonlinearities and large-scale circuits are interesting 
topics to be studied in the future. There is still a problem left. 
In our method, the decoherence will be enhanced if we want 
to generate stronger quantum nonlinearity, which may limit 
our ability to manipulate nonlinear quantum effects, e.g., the 
generation of Schrodinger cat state. This problem is also left 
for the future study. 



Appendix 

Derivation of Eq. HW) : To adiabatically eliminate the de- 
grees of freedom of the squeezed cavity field, we introduce the 
singular perturbation approach ll78l to let k — kq / e, £^ — / e, 
and let e in the final step. In this case, the total system 
composed of the system {S, L, H), the squeezed cavity field, 
and the system {S,Lf,H) can be described by: 



+ y^A{e''t'c + c''e-''t') + ^^/^{L^c 
+ 1 (h^S^Lf - l]Sl) + f (c^S^Lf ^ 



l\sc 



(20) 



From Eq. (30) in Ref. Ii52l and Eq. ( i20] i, the dynamical 
equation of the evolution operator Ut (e) of the total system 
can be expressed as: 



dUt (e) = \{S^ - I) dA 




dt } Ut (e) , 



where 

At- 



/ bl (r) 6i„ (r) dT, Bt^ f (r) dr. 
Jo Jo 



In order to eliminate the degrees of freedom of the cavity mode 
in the singular limit e — ^ 0+, we change to the interaction 
picture by introducing the time evolution operator Vt, which 
satisfies that: 



dVt (e) 




(21) 



Then, we want to consider the evolution of the unitary operator 
Of (e) = Vt (e) Ut (e), which satisfies that 



dUt je) 
dt 



-iH (e) Ut (e) 



(22) 
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The Hamiltonian H (e) is defined by H (e) 

(e) H (e) Vt (e), where 



Hie) = 



-H + -{L'^S'^Lf - VjSL 



Then, we introduce the normaUzed cavity mode in the inter- 
action picture: 



'^vHe)cV,ie). 



From Eq. ( 1211 1. it can be verified that 

dbtie) = 



' 2e 



{L + S'^Lf) dt + SdAt 



(23) 



It can be solved from Eq. ( |23] l that 



h (e) 



[Gi, [t) (c+ct)+G2e (t) (c-ct)] 



+ 



'*o — Co Jo 
+5'dA^ + h.c] 



+ 



+ Co Jo 
+SdAr - h.c] , 



G2e - r) 



- (L + S^Lj) dr 



-{L + S^Lf) dr 



(24) 



where 



Gle (t) = CXp 



G2e it) 



4e 
Ha + Co 



4e 



■ exp 



(^0 - Co) |t| 
2e 

(^o + Co)k| 

2e 



From Eq. (l22l i. we have 



C/t(e) = -i{ [-1(^^-4-5) 



d_^~. 
di 



(25) 

Substituting Eq. (|24] | into Eq. ( l25T l and letting e — > (it is a 
weak convergence in the meaning of Ref. I.78J), we have 



+bl (e) [i {L - S^Lf) - 2Ae~"t'] 
+H+1- (L^S^Lf - l\SL\ \ Ut (e) dt. 



dUt 



- l].S^ (cosh (ro) (L + S^Lf) 



sinhlLt + L^S 



dt 



--Ae'^ ((cosh (ro) + 1) (i + 5+2. 



sinh(ro)(i^ + Lt5)) +h 



dt + i/rft 



1 (Lt - L].S) dBl + ^ (i - ^ti;) dfi 



- (L - 5tL;) rfS/^t ^ 



(26) 



where 

dSJ' = cosh (ro) SdBt + sinh (ro) dSjS't. 
Here, we have used the limit 



Ge(r) = -cxp - — 



5{t), e^O. 



where tr^ (•) is 



Since p = tiE \^(UtPoU}) / ti (UtpaUj^ 
the partial trace over the Hilbert space or the noise dBt, we 
can obtain the master equation (fTOl i of p from Eq. (|26] | with 
simple calculations. Note that the map 

T{t) : po ^ tiE [{UtPoU^) /tr {UtPoU^ 

has good properties, because the final master equation dTOt is 
of the traditional form of squeezing environment which has 
been widely studied in quantum optics ||96l . ||97l . 

Derivation ofEq. dTTl ): Near the optimal point Ug = 1/2, the 
two energy levels of the charge qubit corresponding to n = 0, 1 
are close to each other and far separated from higher-energy 
levels. Thus, in this case, the charge qubit can be looked as a 
two-level system Il85l - ll87l . The Hamiltonian H^t in Eq. JTSl l 
can be written as: 



HgT - -2Ec (1 - 2ng) a, - E^j cos [^-^j 
where the Pauli operators and dx are defined by: 



|0)(0|-|1>(1|, 
|0)(1| + |1)(0|, 



and |0) and |1) are the charge states with the Cooper pairs 
numbers n = 0, 1. 

When Ug — 1/2, the charge qubit is in the charge degenerate 
point. In this case, HgT can be rewritten using the eigenstates 
of the charge qubit as: 



H. 



E°j cos + 77T (a + a^) 

9o 



bl 



where 



and 



-)(- 



->(- 



Since $ 

HgT 



$0/2, we have 



= -E" 



sin -^?7T (a + a^) cos -^7;i„ ( &i„ + 



+ COS ^?7t (a + a^) sin ^»7in (&in + bl^ 
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Expanding HqT to the first order of b\n + b\^^ and the second 
order of {a + a^), we can rewrite Hqx as: 



qT 



= (jJaCi^a — Ej 



Here, we have omitted the term 



which just leads to additional dephasing effects of the charge 
qubit. Assume that the charge qubit always stays in the ground 
state and omit the fast oscillating terms a^^, in the rotating 
wave approximation, we can obtain the effective Hamiltonian 
of the TLR: 



The linear term 



w„ a ' a — 



ttEjTit 



$0 



ttE"", 



$0 

can be compensated by a classical driving field on the TLR, 
and thus the effective Hamiltonian Ht given in Eq. ( fTTI l can 
be obtained. 
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